Optical absorbtion by atomically doped carbon nanotubes 
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We analyze optical absorption by atomically doped carbon nanotubes with a special focus on the 
frequency range close to the atomic transition frequency. We derive the optical absorbtion line-shape 
function and, having analyzed particular achiral nanotubes of different diameters, predict the effect 
of absorbtion line splitting due to strong atom-vacuum-field coupling in small-diameter nanotubes. 
We expect this effect to stimulate relevant experimental efforts and thus to open a path to new 
device applications of atomically doped carbon nanotubes in modern nanotechnologies. 
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PACS numbers: 78.40.Ri, 73.22.-f, 73.63.Fg, 78.67.Ch 

Carbon nanotubes (CNs) are graphene sheets rolled-up 
into cylinders of approximately 1 nm in diameter. Exten- 
sive work carried out worldwide in recent years has re- 
vealed intriguing physical properties of these novel molec- 
ular scale wires 0, y] ■ Nanotubes have been shown to be 
useful for miniaturized electronic, mechanical, electrome- 
chanical, chemical and scanning probe devices and as 
materials for macroscopic composites 0, 0| • Important is 
that their intrinsic properties may be substantially mod- 
ified in a controllable way by doping with extrinsic impu- 
rity atoms, molecules and compounds H, IH @> ■ Re- 
cent successful experiments on encapsulation of single 
atoms into single-wall CNs |8j and their intercalation 
into single- wall CN bundles 0,0], a l° n g with numerous 
studies of monoatomic gas absorbtion by the CN bun- 
dles (see ^3 f° r a review), stimulate an in-depth the- 
oretical analysis of near-field electrodynamical proper- 
ties of atomically doped CNs systems. This is of both 
fundamental and applied interest as it sheds light on 
the peculiarities of the atom-electromagnetic-field inter- 
actions in quasi-lD dispersive and absorbing surround- 
ings, and, eventually, would open routes for new chal- 
lenging nanophotonics applications of atomically doped 
CN systems as various sources of coherent light emitted 
by dopant atoms. 

Typically, there may be two qualitatively different 
regimes of the interaction of an atomic state with a vac- 
uum electromagnetic field in the vicinity of a CN 0, 
HiL IlH llij . They are the weak coupling regime and 
the strong coupling regime. The former is character- 
ized by the monotonic exponential decay dynamics of 
the upper atomic state with the decay rate altered com- 
pared with the free-space value. The latter is, in con- 
trast, strictly non-exponential and is characterized by 
reversible Rabi oscillations where the energy of the ini- 
tially excited atom is periodically exchanged between the 
atom and the field. We have recently shown that the rel- 
ative density of photonic states (DOS) near the CN ef- 
fectively increases due to the presence of additional sur- 
face photonic states coupled with CN electronic quasi- 
particle excitations This causes the atom- vacuum- 
field coupling constant, which is determined by the lo- 



cal (distance-dependent) photonic DOS, to be very sen- 
sitive to the atom-CN-surface distance. At small enough 
distances , th e system may exhibit strong atom-field cou- 
pling giving rise to rearrangement ("dressing") of 
atomic levels by the vacuum-field interaction with forma- 
tion of atomic quasi-lD cavity polaritons. The problem 
of stability of quasi- ID polaritonic states in atomic ally 
doped CNs has been partly addressed in Refs. 0, Q] 
by investigating the atom-nanotube van der Waals in- 
teractions and demonstrating the inapplicability of stan- 
dard weak-coupling-based van der Waals coupling mod- 
els close to the CN surface. In this Letter, we analyze 
optical absorption properties of atomically doped CNs 
with a special focus on the frequency range close to the 
atomic transition frequency and predict the effect of the 
absorbtion line splitting for quasi-lD atomic polaritons 
in small-diameter CNs. 

We consider a two-level atom positioned at the point 
ya near an infinitely long achiral single-wall CN. The 
atom interacts with a quantum electromagnetic field via 
an electric dipole transition of frequency loa- The atomic 
dipole moment, d = d z e z , is assumed to be directed along 
the CN axis (assigned by the unit vector e z ) which is 
chosen to be the ^-quantization axis of the system. The 
contribution of the transverse dipole moment orientation 
is suppressed because of the strong depolarization of the 
transverse field in an isolated CN (the so-called dipole 
antenna effect 0])- The total secondly quantized Hamil- 
tonian of the system is given by |l3|, [lj] 



H = 



dtohuj I dR/ t (R,w)/(R,w) 4 

J dujJdH [gW(r A ,R,w) CT f 
g<->(r A ,R,w)£]/(R,w)+h.c., 



hu>A 



(1) 



with the three items representing the electromagnetic 
field (modified by the presence of the CN), the two-level 
atom and their interaction, respectively. The operators 
f'(R,oj) and /(R, w) are the scalar bosonic field opera- 
tors defined on the surface of the CN (R is the radius- 
vector of an arbitrary point of the CN surface). They play 
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the role of the fundamental dynamical variables of the 
field subsystem and satisfy the standard bosonic commu- 
tation relations. The Pauli operators, & z = \u)(u\ — \1}{1\, 



a = \l) (u\ and a 



t = 



i){l\, describe the atomic sub- 



system and electric dipole transitions between the two 
atomic states, upper \u) and lower \l), separated by 
the transition frequency lua- This (bare) frequency is 
modified by the diamagnetic (^A 2 ) atom- field interac- 
tion yielding the new renormalized transition frequency 
Co A = u A [l - 2/(hLu A ) 2 f™dwfdR\g ± {r Al R,u})\ 2 ] in the 
second term of the Hamiltonian. The dipole atom-field 
interaction matrix elements (ta, R, oj) are given by 
g( ± )=g i ±(w/w j4 )gl where g^H^r^, R, to) = -i{Auj A /c 2 ) 
^TrhujRecj zz {tu) d z ±m G zz {r A ,R,oj) with X ^G ZZ be- 
ing the zz-component of the transverse (longitudinal) 
Green tensor (with respect to the first variable) of 
the electromagnetic subsystem and <t zz (uj) representing 
the CN surface axial conductivity. The matrix elements 
g-L(||) nave the property of JdR \g ± ^(rA, R, uj)\ 2 = 

(H 2 /2Tr)(LUA/L0) 2 r a (L0)^(T Al L0) with ^(t A> u) = 

Im ±( -"- ) G' z l z '"' ) (r J 4,r J 4,a;)/ImG° z (a;) being the transverse 
(longitudinal) distance-dependent (local) photonic DOS 
functions and Tq(lli) — 8Truj 2 d 2 ImG® z (uj)/3frc 2 represent- 
ing the atomic spontaneous decay rate in vacuum where 
In\ G° zz (lu) —lu/6ttc is the vacuum imaginary Green ten- 
sor zz-component. The Hamiltonian involves only two 
standard approximations: the electric dipole approxima- 
tion and two-level approximation. The rotating wave ap- 
proximation commonly used is not applied, and the dia- 
magnetic term of the atom-field interaction is not ne- 
glected (as opposed to, e.g., Refs. 0,0,0). Quantum 
electrodynamics (QED) of the two-level atom close to the 
CN is thus described in terms of only two intrinsic char- 
acteristics of the electromagnetic subsystem — the trans- 
verse and longitudinal local photonic DOS functions. 

In deriving the absorbtion/emission spectral line sha- 
pe, we follow the general quantum approach of Ref. |l7j |. 
(Obviously, the absorbtion line shape coincides with the 
emission line shape if the monochromatic incident light 
beam is being used in the absorbtion experiment.) When 
the atom is initially in the upper state and the field sub- 
system is in vacuum, the time-dependent wave function 
of the whole system can be written as 

\m) = aWe-^^'MKO}) (2) 

dR / du;Ci{R,uj,t)e- l{u - CjA/2)t \l}\l(R,oj)), 



where |{0}) is the vacuum state of the field subsystem, 
|{1(R, u>)}) is its excited state with the field being in the 
single-quantum Fock state, and C u and Ci are the popu- 
lation probability amplitudes of the upper state and the 
lower state of the whole system, respectively. Applying 
the Hamiltonian to the wave function J2J), one has 



C„(t) 



dRg< + >(rx,R,w) 



(3) 



Q{R,LU,t): 



Ci(R,oj,t)e 



g(+\r A ,R,u>)YC u (t)e 



-i{uj~Cj A)t 



i(uj—uJA)t 



(4) 



In terms of the probability amplitudes above, the emis- 
sion/ absorbtion spectral line shape I(lu) is given by 
JdR \Ci(R, uj, t — ► oo)| 2 , yielding, after the substitution 
of Ci obtained by integration [under initial conditions 
Q(R,a;,0) = 0, C„(0) = 1] of Eq. ©, the result 
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I(x)=I (x) 



dr C„(t) e 



i{x — x 



(5) 



where I (x) = {f (x)/2Tr}{(x A /x) 2 £- L (x) + ^(x)} and we 
have, for convenience, introduced the dimensionless vari- 
ables / = 2j I/h, T = ?iro/27o, x = huj/2'f , and 
t = 2-f t/h with 70 = 2.7 eV being the carbon nearest 
neighbor hopping integral appearing in the CN surface 
axial conductivity a zz . 

The upper state population probability amplitude 
C u (t) in Eq. (JSJ) is given by the Volterra integral equation 
[obtained by substituting the result of the formal integra- 
tion of Eq. (@J into Eq. with the kernel determined 
by the local photonic DOS functions ^\r A , x) [HE1- 
Thus, the numerical solution is only possible for the line 
shape I(x), strictly speaking. This, however, offers very 
little physical insight into the problem of optical absorb- 
tion by atomically doped CNs under different atom-field 
coupling regimes. Therefore, we choose a simple analyti- 
cal approach valid for those atomic transition frequen- 
cies xa which are located in the vicinity of the reso- 
nance frequencies x r of the DOS functions Within 
this approach, ^- L CH) (r^4 , x~x r ) are approximated by the 
Lorentzians of the same half- width- at-half-maxima Sx r , 
thus making it possible to solve the integral equation for 
C u analytically to obtain 0, 0] 



1 



Sx r 



C u (r)w- 1 + - | e 

u{ ' 2\ ^5x 2 -X 2 

1 / 8x r 



2 \ ^Sx 2 - X 2 



-^(SXr + ^/Sxl-X^T 



with 



X 



\J 2Sx r r (xA)^- L (r A , xa) ■ 



(7) 



This solution is valid for 53^4 ~ x r whatever the atom-field 
coupling strength is, yielding the exponential decay of the 
upper state population, |C u (r)| 2 w exp[— T(xa)t] with 
the rate T Tq^, in the weak coupling regime where 
(X/8x r ) 2 <C 1, and the decay via damped Rabi oscilla- 
tions, |C„(t)| 2 m exp(— 5x r r) cos 2 (Xr/2), in the strong 
coupling regime where (X/Sx r ) 2 1. 

Substituting Eq. @ into Eq. © and integrating 
over t, one arrives at the expression 



I(x) = I Q (x A ) 



(x ~ xa) 2 + 8x 2 . 



[(x - x A ) 2 - X 2 /4} 2 + 6x 2 (x - x A ) 2 



(8) 
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representing the emission/absorbtion spectral line shape 
for frequencies x ~ x r ~ xa regardless of the atom-field 
coupling strength. The line shape is clearly seen to be of 
a symmetric two-peak structure in the strong coupling 
regime where (X/5x r ) 2 ^> 1. The exact peak positions 

are x lj2 = x A ± (X/2)J '^/T+ 8{Sx r /X) 2 - 4(Sx r /X) 2 , 
separated from each other by x\ — X2 ~ X with the 
Rabi frequency X given by Eq. @. In the weak cou- 
pling regime (X/8x r ) 2 <C 1, and x\,2 become complex, 
indicating that there are no longer peaks at these fre- 
quencies. As this takes place, Eq. (jHJ is approximated 
with the weak coupling condition, the fact that x ~ xa 
and Eq. J7J, to give the well-known Lorentzian I(x) = 
Io(x A )/[{x- x A ) 2 + f 2 (x A )/4] of the half-width-at-half- 
maximum T(xa)/2, peaked at x = xa- 

To compute the absorbtion line shapes of particular 
CNs, one needs to know their transverse local photonic 
DOS functions ^(r^,^) in Eqs. and ijHJ. These are 
basically determined by the CN surface axial conductiv- 
ity <Jzz- The latter one was calculated in the relaxation- 
time approximation (relaxation time 3 x 10 -12 s) at tem- 
perature 300 K. The function was calculated there- 
after as described in Refs. 0,0,EI0]- The vacuum 
spontaneous decay rate was estimated from the expres- 
sion Tq(x) — a 3 x (a — 1/137 is the fine-structure con- 
stant) valid for atomic systems with Coulomb interac- 
tion 0. 

Figure^shows the computed £ (x) and normalized op- 
tical absorption spectra A(x) = I(x)/ I pea k for the atom 
close to the (9,0) CN. When the atom, being outside, ap- 
proaches the CN surface, the function £ goes up [shown 
in Fig. 1(a); the vertical dashed lines indicate the DOS 
peak frequencies for which A(x) in Figs.^b) and^c) are 
calculated] due to the increasing role of additional surface 
photonic modes of the nanotube. This causes the atom- 
vacuum-field coupling to increase with the effect [shown 
in Fig. 1(b)] of first broadening and then, when the strong 
atom-field coupling regime is realized, splitting the ini- 
tial Lorentzian line into two symmetric components sep- 
arated by the Rabi frequency 1x0.1 corresponding to 
the energy O.lx 270 = 0.54 eV. In Fig.QJc) much stronger 
Rabi splitting of X ss 0.3 is demonstrated (Rabi energy 
0.3 x 270 = 1.62 eV) for the atom located in the center 
of the CN. This is in agreement with the fact that, due 
to the nanotube curvature, the effective interaction area 
between the atom and the CN surface is larger when the 
atom is inside rather when it is outside the CN |l3l fl4| . 

Figure l^a) presents ^(x) for the atom in the center 
of the (5,5), (10,10) and (23,0) CNs. It is seen to de- 
crease with increasing the CN radius, representing the 
decrease of the atom-field coupling strength as the atom 
moves away from the CN wall. To calculate the normal- 
ized absorbtion curves A(x) — I (x) / I pea k in this case, we 
have fixed xa = 0.45 [indicated by the vertical dashed line 
in Fig. Ola)] because this is the approximate peak posi- 
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FIG. 1: (Color online) Transverse local photonic DOS (a) and 
normalized absorbtion line shapes (b) for the atom at different 
distances outside the (9,0) CN. (c) Normalized absorbtion line 
shape for the atom in the center of the (9,0) CN. 



tion of £ for all the three CNs. The result is shown in 
Fig. [2b). For the small radius (5,5) CN, the absorbtion 
line is split into two components, indicating the strong 
atom-field coupling with the Rabi splitting X«0.2 cor- 
responding to the energy of 0.2 x 270 = 1.08 eV. The 
splitting is disappearing with increasing CN radius, and 
the absorbtion line shape becomes Lorentzian with the 
narrower widths for the larger radii CNs. 

Recently, Rabi splitting ^140 — 400 /xeV was detected 
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FIG. 2: (Color online) Transverse local photonic DOS (a) 
and normalized absorbtion line shapes (b) for the atom in the 
center of the CNs of increasing radii. 



in the photoluminescence experiments for quasi-OD exci- 
tonic polaritons in quantum dots in semiconductor mi- 
crocavities [lji |2fJ, |2J ■ The effect we report in this Let- 
ter is at least 3 orders of magnitude larger. This comes 
from the fact that the typical atomic binding energies, 
are at least 3 orders of magnitude larger than the typ- 
ical excitonic binding energies in solids. In terms of the 
cavity-QED (see, e.g., Ref. 22\), the coupling constant 
of the atom-cavity interaction is given in our notations 
by hg = d z (2nhujA/V) 1 ^ 2 with V being the effective vol- 
ume of the field mode the atom interacts with. In our 
case V~ i:R 2 cn (\ A /2), which is — 10 2 nm 3 = 10~ 7 /im 3 for 
the CNs with diameters ~ 1 nm in the optical spectral 
region (as apposed to ~ 10~ 1 -^10~ 2 ton 3 for semiconduc- 
tor microcavities in Refs. [HJ, |2fJ, |2lJ ) . Small mode vol- 
umes yield large atom-cavity coupling constants in small- 
diameter CNs. For example, in the situation shown in 
Fig. |2Ja) one has %~0.3 eV for the atom in the center 
of the (5,5) CN [R cn = 0.34 nm, d z ~ er ~ e(e 2 /hu> A )], 
whereas the "cavity" linewidth, which is related to the 
Purcell factor ^(u) a) via the mode quality factor when 
ujA~uj r H3, is h'y c = 6irhc 3 /u 2 a ^(£j a )V ~0.03, so that 
the strong coupling condition g/"f c 3> 1 has been well 
satisfied. 



To summarize, we have shown that, similar to semicon- 
ductor microcavities and photonic band-gap materials, 
carbon nanotubes may qualitatively change the charac- 
ter of the atom-electromagnetic- field interaction, yielding 
strong atom-field coupling. We predict the effect of the 
optical absorbtion line splitting of ~ 1 eV in the vicinity 
of the energy of the atomic transition for quasi-lD atomic 
polaritons in small-diameter atomically doped CNs. We 
expect this effect to stimulate relevant experimental ef- 
forts and thereby to open a path to novel device appli- 
cations of atomically doped CNs. For example, quasi-OD 
excitonic polaritons in quantum dots in semiconductor 
microcavities were suggested recently to be a possible 
way to produce the excitonic qubit entanglement [23| . 
Leaving the details for a forthcoming publication, we con- 
clude that, being strongly coupled to the (resonator- like) 
cylindrical nanotube environment, the two atomic quasi- 
lD polaritons can be entangled as well, thus challenging 
novel applications of atomically doped CNs in quantum 
communication technologies. 
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No W911NF-05-1-0502 and NAG3-804. 
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